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1 Introduction 

We investigate a system being time discretization of two dimensional Navier-Stokes 
equations in the isentropic regime 

^ (g k - g^ 1 ) + div(g k v k ) = 
_i_ (gk v k _ gk-i v k-i^ _j_ div(g k v k (g> v k ) - fiAv k - (/i + z/)V divv k + Vn(g k ) = 0, 

where C f 2 is a fixed domain, v k : — > M 2 - the velocity field, g k : f2 — > Wq- the 
density, it : Mjj" — > M- the internal pressure given by the constitutive equation 

*{Q k ) = (Q k T, 7 > 2. 
We assume that the walls of Q are rigid and that the fluid slips at the boundary 



v k - n = 0, at dtt 
n-T(v k ,TT)-T + fv k -r = at dtt, 



(2) 
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where T(v k ,ir) = 2fiD(v k ) + [y divu — tt)I. 

The conditions ([2]) are known as the Navier or friction relations which means, that 
unlike in the case of complete slip of the fluid against the boundary, the friction effects, 
described by / > 0, may also be present. The customary zero Dirichlet condition may 
be understood as a special case of the above, when / — > oo. 
We will always assume that our initial conditions g°,v° satisfy 

a.e. %n Q, g° G L 7 (ft), 
A°eI 27/(7+1) (Q), ^°(f ) 2 g Li(Q). 

The first main goal of this paper is to show that for At = const, the solutions of such 
a system exist in a sense of the following definition. 

Definition 1. We say, the pair of functions (g k ,v k ) G L 7 (fl) x W^ity, v k - n = at 
dQ is a weak solution to (TJp-([^) provided 

J g k v k -V v dx = ^- t J{g k - g k ~ 1 )(p dx, G C°°{Q), 



n 



and 



2^ J (g k v k - g k - l v k - l )if dx- J g k v k ® v k : Vtp dx + 2fi J T>(v k ) : D(^) dx 

n n n 

+ v J divv k divip dx - J n(g k ) divtp dx + J f(v k ■ r)((p • r) dS = 0, 
n n an 

Vip G C 00 ^); cp-n = at dn. 

The first main result reads as follows. 

Theorem 1. Let Q G C 2 be a bounded domain, At = const., \x > 0, 2/x + 3u > 0, 
7 > 2, / > ; g k ~ 1 > 0. Then there exists a weak solution to (QP-([J|) such that 

g k G Loo(fi) and g k > 0, 
v k G M^(n) Vp < oo, 

/(ix = / g k ~ l dx, 



moreover ||^ fc ||oo < (At) 3 ^ 1 

The result we present here was already stated without requiring any assumption 
on the smallness of initial data f.i. in the monograph of Lions [5] for the zero Dirichlet 
condition when Q is bounded and for the whole space. It was used there as a tool in 
analysing the steady and non-steady cases. The approach presented there or in other 
works based on the Feireisl idea [S], [2], [I] benefit from the properties of the effective 
viscous flux. Our technique allows for essential reduction of the number of technical 
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tricks and enables to get the required regularity for the density directly at the level 
of approximate system. 

The second result refers to a passage to the limit with length of time interval At — > 0. 
We will show that for such a case our solution tends to the weak solution of non-steady 
compressible Navier-Stokes system with a slip boundary condition: 

Qt + div(fw) = in ft 
(gv)t + div(£w <g> v) — [iAv — (p+ u)V divv + Vn(g) = in Q , , 

v . n = o at on ^ 

n- T(v, it)- t + fv- t = at dfl, 

in sense of the following definition. 

Definition 2. We say, the pair of functions (g,v) G Loo(L 7 ) x /^(W 7 ^), v ■ n — at 
dQ is a weak solution to (TJ|) provided 

T 

[gift + gv ■ V<p) dxdt = 0, G C 6 °°([0, T)xfi), 




o n 

and 

T 

(gvipt + gv <S> v : V x ip + ir(g) div x <^) dxdt = 

o n 

T T 





(2/iD x (t>) : ~D x (ip) + v div x t> div^) dxdt + J j f(v ■ r)(tp ■ r) dSdt, 
on o an 

Vy? G C c °°([0, T) x fi); 9? ■ n = at (4) 

The existence of solutions to the non-steady system is provided by our second main 
result. 

Theorem 2. Under the hypotheses of Theorem^ the solution (g k ,v k ) converges to 
(g,v) as At — > + weakly (weakly*) in L^L^) x L 2 (H / 2 1 ). 

Moreover g belongs to L 1+ i(Q x (0, T)) and the following energy inequality is satisfied 
for almost all t G [0, T] 



Of 2 



(T)dx + ^-jJ g f (T)dx + J J (2{i\D(v)\ 2 + v( divv) 2 ) dxdt 
n on 

T 



+ J J f(vT) 2 dx dt < C(g u ,v°). 
o an 
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In the following section we will show the existence and uniqueness of regular solution 
to the problem being the new e— approximation scheme for the time-discretized Navier- 
Stokes equations. Although the proof is based on the standard fixed-point method, we 
will precisely present most of steps in view of the fact that our approximation affects 
the nonlinear term too. Our solution (g k , v k ) will be obtained as a weak limit as 
e — > + of the sequences (g k , v k ). This limit process will be carried out in Section 3 
by using some uniform estimates and the following property of the density sequence 

lim {i 6 O : g k (x) > m}| =0 
for m sufficiently large, which enables to show the convergence of the pressure. 



2 Approximation 

In this section we present a scheme of approximation being a modification of the one 
introduced by Mucha Pokorny |6J for the steady case. It is needed to investigate the 
issue of existence of solutions in the case when the time step (At) is cosnstant and while 
disposing a sufficient information for the density and velocity at the k — 1 moment of 
time. Although for further purposes there is a necessity to keep trace of the dependence 
on these quantities in almost all estimates. 
Let 

a = J - 

Ai' ftj\ 

h = £ fc ~\ g = g k , v = v k , g = v k ~ l . 

The objective of this part of work will be then to examine the following approximative 
system: 

a(g- hK(g)) + div(K(g)gv) - eAg = 
a (gv - hg) + div(K(g)gv <g> v) - [iAv - (n + v)V divw + VP(g) + eV^Vf = 

§£ = at dQ, (6) 
v ■ n — at dQ, 
n-T(v,P(g))-r + fv-r = at dQ, 

we will write simply g, v istead of g € , v t when no confusion can arise. The other deno- 
tations are the following: 

e 



P(g) = 7 J s ^ l K{s)ds, 







g < mi, 

K(g) = { g > m 2 , 




where 



g G (m 1 ,m 2 ), 
and 

K(-) E C 1 (1R) K'(g) <0 in ( mi ,m 2 ), 

for some constants m 1; m 2 . To avoid the difficulties conected with the case when 
mi m2 we set the difference 777,2 — Wi to be constant, equal 1. 
The existence of a regular solution is guaranteed by the theorem. 
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Theorem 3. Let Q G C 2 , e, go, > 0. Then there exist a regular solution (g,v) to 

(EP, q G % 2 (^), u G W^ 2 (fi) /or a//p < oo. 

Moreover 

< q < m 2 m fi, (7) 
^dx < / /icfe. (8) 



Proof. We assume, that £>, t> are regular solutions to ([2]) and prove some estimates first, 
after we go on with the existence. 
Step 1. Proof of (TJj. 

Integrating the first equation of ([6]) over Q one gets 

a j ' (g — hK(g))dx + / K(q)qv ■ ndS — e J ~^~dS = 0, 
n an an 

the boudary integrals vanish and due to the definition of K(-) we truly have 

gdx = / K(g)hdx < / hdx. 



n n n 

Step 2. Non-negativity of g. 

Assume, that we have h > in Q, the proof follows by the induction. We integrate 
first equation of ([6]) over Q~ = {x G Q : g(x) < 0} 

a J (g — K(g)h)dx+ J K(g)gv-ndS — e J -^-dS = 0, 

n- 8U- au- 

the first boundary integral vanishes since either g or v ■ n equals at dQ~ . Moreover, 
we know that > at dQ~~ , hence 



but this leads to conclusion that = and consequently £ > in fl. 
Step 3. Upper bound for g. 

Assume that h < m2. This time we integrate the approximate continuity equation over 
£}+ = { x G Q : £>(a;) > m 2 } 

a J (g — K(g)h)dx + J K(g)gv-ndS — e J ^-dS = 0, 

n+ an+ an+ 

At dn + we have fj[ < and either iT(^) or f • n equals 0. Thus, in the similar way as 
previously, the observation 



J gdx < m 2 J K{g)dx < 



n+ n+ 
5 



implies that g < m 2 in Vt. 
Step 4. Existence. 

In accordance with our denotations the proof of existence of approximate solutions is 
almost identical to the one presented in [6]. In the first step we define for p £ [1, oo]: 

M p = {w E Wj(fi); w ■ n = at dtt} . 

and we claim that the following proposition, which is the analogue of Proposition 3.1. 
from [6 J holds true. 

Proposition 4. Let assumptions of theorem [3] be satisfied. Then the operator S : 
Moo -> W%{&), where 

S(v) = g, 

ag + div(K(g)gv) — eAg = ahK(g) in Q 

|^ = at on 

is well defined for any p < oo. Moreover 

• Q — S(v) satisfy 

S edx - 

n 

• If h > i/ien g > a.e. in fL 

• ^/ |l v l|i,oo < L, L > i/ien 

IMkp < C(e,p,n)(l + L)||%, Kp<oo, (9) 

The only difference in the formulation and the proof relates to the fact that h is 
not a constant parameter any more and that there appears g instead of v. But the 
assumtion that the regular solution in k — 1 moment of time exist allows to replace the 
modulus by the L p norm of h. 

In the next step of proof of Theorem |3] we will consider the Lame operator 

T : ^ Moo 

defined as follows: w = T(v) is a solution to the problem 

—/lAw — ((/, + u)V divw = ahg — agv — div(K(g)gv (8) v) — VP(g) — eVgVv = 

= F(g,v,h,g) 
w ■ n = at dQ, 
n- (2p,D(w) + v divw/)- r + fv ■ r = at Oil 

(10) 

Employing the Larey-Schauder fixed point theorem for the operator T we can almost 
rewrite the proof of analogous fact [8] or [6]. The only part that that deserves more 
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careful study is the energy estimate. This in turn together with some information 
about the pressure P(g) will enable to pass to the limit with the length of time interval 
At. First observe that ffTOjl i with w — v and g = S(v) holds with a solution itself as a 
test function, therefore 

a J gv 2 + J dxv(K(g)gv <g) v)v — fj, J (Av)v — (/x + v) J (V drvv)v + JvP(g)v 
n n n n a 



+ e J VgVvv = a / hgv. 
h 

Next, integrating by parts and using condition on the boundary 
a J Q v2 + \J div (K(g)gv)v 2 + 2 fi J \D(v)\ 2 + u J div 2 v + J f(v-r) 2 



n n an 

— f div(K(g)gv)g' y - 1 - | J Agv 2 = a J hgv, 



7-1 



including the information contained in (|6j) i one gets 
-a f(g + K(g)h)v 2 + 2fi [ \D{v)\ 2 + v [ div 2 v + I f(v-r) 2 



2 



<-m 



hgv, 



+ — jya J g 1 ~r~[ a J Q 1 ' 1 K(g)h + 1 t j \Wg\ 2 g^ 2 = a J, 

a n a a 

now we add end substract \ot j n hg 2 

J {gv 2 - hg 2 ) + -a J h\v - g\ 2 + 2/x ^ |£>0)| 2 + f ^ div 2 v 
n n n n 

+ y /(^) 2 + 7^7« / ^-^37" / £? 7 ~^(^ + 7e y |V^| 2 <0, (11) 



an n 
next we add and substract — f n ft, 7 

7—1 JS2 



-a y (gr; 2 - /u/ 2 ) + it* y %-s| 2 + 2/x y \D(v)\ 2 + u j div 2 v + j f(vr) 2 
n n n n on 

+ ^Y« y (^ 7 -^ 7 ) + ^Ta y ((7-l)^ 7 + ^ 7 -7^ 7_1 ^)^)+7e J |V^| 2 < 0. 



n 

(121 



Note, that since g, h>0 and < 1 we have that (7 — l)f5 7 + /i 7 — 7£ 7 1 K(g)h > 0. 

Referring to our oryginal denotation we may now sum (jl2l) from = 1 to = n and 
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obtain the following bounds: 



sup T a\\^\\y + -a\\g n (v n )% < T a\\gX + o<x\\e (v°) 2 \\i, 

0<n<M 7 ~ 1 J I 7 - 1 '2 



thus 



sup WQX+ \\e n (v n ) 2 \\i<C(g ,v°^,n), (13) 

0<n<A/ 



in particular C is independent of k, e and a, moreover 

M 



J2 J [f? fc ~V - v k - l \ 2 + (7 - l)(g k y + ( e k ~ 1 y - ^(Q k ) ry - 1 K( Q k ) e k - 1 ] < C (14) 



k=1 n 



with the same constant C. The information contained here turns aut to be one of the 
crutial importance at the second stage of this work while showing that the passage 
with At — > gives the solution to the evolutionary case. Namely, since for 7 > 2 there 
exists a positive constant 5, such that 

( 7 - i)( g k y + (g^y - i^y^K^g*- 1 > 5\ Q k - Q k ^\\ 

hence (Tl~4"|) ensures 

M 

^\Q k ~ Q k ~ l V <C. (15) 



Additionaly we have 



and by Korn's inequality 



k=l 



M 



J2\\Dv k \\l<aC 



k=l 



M 



Y,\\A\l2<t*C (16) 



k=l 



here the constant C depends also on /1 and v. 
Finally we also get 

M 

£liv(^)^<-a (17) 



k=l 



This information allows us to repeat the procedure described in [8], which together 
with the Proposition H] yield the existence of regular solutions, and hence the proof of 
Theorem [3] is complete. 

Apart from the information resulting from the first a priori estimate , the limit passage 
requires also some estimates independent e, a and m,2- 

First of them is the estimate for the norm of gradient of the density. Observe that 
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multiplying ©1 by £> and integrating over Q one get 



6 / \\7 Q \ 2 = a / hK(g)g-a / g 2 - / K{g)gv-Vg 



n n n n 





< «Cm 2 + y v • V j J K(t)t dt J = aCm 2 - J divv j J K(t)t dt 

fi \0 / fi \0 y 

< «Cm2 + J I divf |^i 2 < aCrri2 + -y/aCm 2 . 
h 

This means that ||Vf?||2 may blow up as e — > + , however we can provide that e||V^||2 
will tend to zero, i.e. 

e\\Vg\\2<V~eC{a,m 2 ), (18) 
for some constant C independent of e. 

Now we would like obtain integrability of the pressure with the power 2, as previously 
independently of e and, if possible, of 777.3. 

Therefore the choise of an appropriate test function seems to be obvious: 

$ = b(p(q)-{p(q)}), mn 

$ = at dQ 

where B is the Bogovskii operator. By Lemma 3.17 from [S] and the Poincare inequality 
we have: 

< c(p,Q)\\P(g)\\ p , ||V$|| P < c(p,Q)\\P(g)\\ p (19) 

(ft ^/P<2 
< p < 00, p = < arbitrary > 1 if p = 2 
[00 if p > 2. 

From this testing, the following identity appears: 
/ \ 2 

,2 1 



P{qY = — - I / P{g) I + a I (gv-hg)$+fi / : V$ + (ji+u) / divw div$ 



n n 



- / K{g)gv <g) u : V$ + e / V^Vw$ = Ij. 
^ ^ i=l 



Now each term will be estimated separately. 

(i) By the estimate (I13p and the definition of P the first one comes strightforward 



(ii) The relation (1191) together with the estimate (ITS]) imply 

I 2 = aj(gv- hg)<$> dx < Ca (h^Vf^ + \\h\\ 1 P\\h?\\ 1 1 /2 ) \\P(g)\\ 2 
n 

< Ca\\P{g)\\ 2 . 

(iii) We also have ||V$||2 < || £? 7 || 2, thus 

h + h = V J VwV$ + (// + v) J divv div$ < C\\v\\ 2 \\P{e)\\2 
n n 

<C a ^ 2 \\P(g)\\ 2 . 

(iv) Since the modulus of K is less than 1, the Holder's inequality and imbedding 
mentioned above lead to 



h = 
n 



J K(g)gv <g> v : V$ < C|M| 7 |M|? >2 ||P(f?)|| 2 < Ca\\P(g)\\ 2 . 



(v) Finally, epmloying the Holder's inequality we may get that 

J e = e / VeV«$<e||Ve|| g ||T;||i, 2 ||P(e)||2, 



n 



for some q > 2. To get the estimate for ||V^|| g we need to interpret the approximate 
continuity equation as a Neumann-boundary problem 

— eA^ = divfr in O 



^ = b-n at 0^, ( 2 °) 
with the right hand side 



b = aB(K(g)h — g) — K(g)gv. 

From the classical theory we know that if dQ is smooth enough and if b G (L p (Q)) 2 , 
then there exists the unique g G Wp(Q) satisfying f l20|) in the weak sence, such that 
j n gdx = const. Moreover 

l|V^<^%|| P . (21) 
Now, in our case assume that 7 > q > 2 then the q-norm of b may be estimated as 

11% < afllell, + \\h\\ q ) + Cy| 7 |Mli,2 < C ia + C 2 ^, (22) 
thus the observation (1211) yields the following 

h = e J VgVv<$> < (da 3/2 + C 2 a)\\P(g)\\ 2 . 
h 
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Gathering the estimates terms ij for i = 1, . . . , 6 one can easily see that 

\\P(g)\\ 2 <Ca^ 2 , (23) 
where the constant C does not depend on e nor 777.2. 

Now our aim will be to estimate the norm of Vv in L q (Q) for q > 2. For this puropse 
we will apply to the system (TT01) the following Lemma (for the proof, see [6] Lemma 
3.3.). 

Lemma 5. Let 1 < p < 00, £1 G C 2 , F G (M 2p /(p + 2))*, /i > 0, 2/i + 3z/ > 0. T/ien 
i/iere exists the unique w G M p , solution to DD\) . Moreover 

\\w\\ hp <C(p,Q)\\F\\ (Mp/(p _ 1)r . 

IfQe C l+2 , F G Wj(fi), / = 0, 1, . . . £/ien to G Wj +1 (^) and 

||«;|| i+2 , p <C(p,^)||F||^. 

If we consider the approximate momentum equation as a part of Lame system with 
w = v we will get the estimate for the norm of Vv in L q (Q) 

\\Vv\\ q < C(a\\gv\\ 2q/{q+2) +011^1125/(5+2)) + \\K(g)gv®v\\ q + \\P(g)\\ q 

+ e||V^Vw|| 2 g/(g + 2)). 

Recalling 7 > 2, by (TT3]) and by (TT6I) one gets 

a\\gv\\ 2q /( q +2) + a\\hg\\ 2q /( q+ 2) < Ca(||^|| 2 + \\hg\\ 2 ) < Ca 3/2 . 
By the definition of P and the Holder's inequality we also have 

\\K{q)0v9v\\ 9 < C\\P(gW q/ Jv\\l 2 < Ca\\P{g)\\ l £. 

At this step there is a need to include the estimates depending on the parameter m 2 , 
more precisely we will use 

\\P(Q)\U < ll^ll^i^ll^ < Ca^mt 2/qh , 
e\\V gVv\\ 2q /( q+2 ) < e||V^|| g ||u||i, 2 < Ca 3/2 m 2 , 

where the last inequality is obtained by replacing in (1221) the norms of g by ||^||oo < m 2 
if q > 7; for q < 7 we can use the estimate for the L 7 -norm of g. 

Summarising, we have shown that ||Vf|| g < C(m 2 ,a) with a constant C(m 2 ,a) inde- 
pendent of e. Particulary for 2 < q < 7 one has 

||Vv||, < C(a + c?l 2 + a 3 /«mt 2/qh )- (24) 
Before passing to the zero limit with e we will compute a priori estimate of the vorticity 

dv 2 dvi 



cu = votv 



dx\ dx 2 
11 



Differentiating n • v = at d£l with respect to the length parameter and combining it 
with the last boundary condition in the system (J2]) we obtain: 

u = \ 2x - v ■ t at dfl. 

Taking the rotation of (jHJta, we get 

— fiAu = —aiot{hg — gv) — rot div(K(g)gv ® v) — erot(V^Vt>). (25) 

Denote 

u) = ui + oj 2 , (26) 

where uj\, u 2 satisfy: 

—fiAoJi = —rot div(K(g)gv <g) v) in Q, 
U\ = at dQ, 

—fiAu2 = —aiot(hg — gv) — erot(V gVv) in Q, 
ui2 = \2x — — J v • r at dVt. 



For the weak solutions coi,u 2 of the above problems one get the following estimates: 

ll^illg < C\\K(g)gv <g> v\\ q < Ca 

where for q < 7, C is independent of m 2 and for q > 7, C = Com 2 

\\u2\\i, q < C(a\\hg\\ q + a\\gv\\ g + e\\VgVv\\ q ) + C(fi)\\v ■ r||i_ 1/9i9j9n , 

thus for q < 2, the Holder's inequality, the imbedding W 2 ^ 2 {dVt) C Wq 1 ^ q (dQ) and 
the trace theorem imply 

||w 2 ||i,g < C(a\\hg\\^ + a\\gv\\^ + e\\Vg\\ 2q /(2- q )\\Vv\\2) + C(tt)\\v\\ h2 

7+1 7+1 

< C(a + a 2 m 2 ) + C(tt)a 1/2 , 
for q > 2 we must use m 2 -dependent estimates of gradient of v in higher norms, thus 

11^2 1| l,^ < C(a,m 2 ) 
and the dependence of m 2 is higher then linear. 

3 Passage to the limit 

This section is devoted to the passage with e — > in the system Recall that so far 
we have obtained the following estimates: 

||£ e ||oo<w 2 , |K||i,2 < Ca, (27) 

\\P(Qe)h<Ca 3 / 2 (28) 

lkl|i l(? + e 1/2 ||V£ £ || 2 < C(m 2 ,a) for 1 < q < 00, (29) 

e|| Vg t Vv t \\ q < C(m 2 , a) for 1 <q < 00. (30) 
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The two last estimates together with the interpolation inequality imply that for 5 
sufficiently small we additionally have: 

e 1_<5 ||V^Vv e ||g < C(m 2 ,a) for 1 < q < oo. 

Therefore, at least for an appropriately chosen subsequence: 

g e ^* g in L^Q), 

P(Q.)^PU) mL 2 (ft), 

v e ^v in W^ft), 

e\/g € ->■ in L 2 (fi), 

eV g e Vv e — > in L q (Q) for 1 < q < oo, 

where the line over a term denotes its weak limit. 

These information allow us to pass to the limit in our approximative system: 



a 



g - hK(gf) + div{K{g)gv) = 



a (gv — hg) + div(K(g)gv (8) v) — /j,Av — (// + ^)V divt> + VP(g) = 

v ■ n — at <9f2, 
n-T(u,P(g))-r + /u-r = at 0ft. 

To show that we have realy found the solution to our initial problem there left several 
questions that need to find the answer. 

Firstly, if we can get rid of K(g) that remains at several places, i.e. if we can prove 
that K(g) = 1 a.e. in fl This, as we shall see below, is equivalent with showing that 
there can be suitably chosen constant m sufficiently large but still sharply smaller than 
the it a priori bound for a density, such that the measure of the set 

{x E ft : g en (x) > m} 

tends to zero for some subsequence e n — > + . Indeed, as this implies that for any 
smooth function n one get 

J Qe n K(g en )n dx = J g tn n dx + J (K(g e J - l)g tn n dx. 

n n {ft„>mi} 

If we choose m\ sufficiently close to m 2 and additionally assure that m <m x then the 
last term on the right hand side disappears as e n goes to 0, and thus we truly have 

hm / Qe n K(g en )n dx = [ gn dx, \/n e C°°(ft). 

fn^0+ J J 

n n 

The next difficulty concerns the convergence in the nonlinear term i.e. is it true that 
P(g) = P{g)- The positive answer can be obtained in a rather standard way, and 
at the stage when one already knows that K(g) = 1 it reduces to proving the strong 
convergence for the density sequence. 
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Finally, what does the condition @H])4 mean, in other words in which sense is it sat- 
isfied? Having solved the two previous problem this is quite easy to see that this 
boundary condition can be recovered while passing to the limit in a weak formulation 
corresponding to the momentum equation. 

Now our aim will be to precisely justify the considerations developed above. For this 
purpose we will adapt a kind of technique widely used for these type of problems, more 
precisely we will take advantage of some properties of the effective viscous flux denoted 
in this paper by G. 

Introducing the Helmholtz decomposition of the velocity vector field defined as: 

v = V0 + V ± A, (32) 

where the divergence-free part V ± A = gf^-J A and the gradient 

by: 

J AA = rotf in Vt J A(f> = divv in Q, 

\ V^A = at dn ' \ = at dtt ' 

we can transform the limit equation (IBTjl o into the form: 

VG = ahg - agv - div(K(g)gv ® v) + /iAVM, (34) 

where VG = V f-(2// + v) A<p + P{g)j . By the observation f n Gdx = J n P(g)dx < 
oo, we control the mean value of G and thus the expression 

G= (2n + u)A4> + Pjg) 

may be accepted as a correct definition of G. 

Due to (133|) . and the classical theory for the laplacian supplemented by the Neuman- 
boundary condition 

\\G\\ 2 <C(\\Vv\\ 2 +\\Pjoj\\ 2 ) <C(a). 

The next goal is to show the boudedness of the norm of G. By the fact that the 
mean value of G is controlled we can employ the Poincare's inequality and the Sobolev 
embedding theorem it is sufficient to prove the following fact: 

Lemma 6. For q > 2 we have: 

\\VG\\ q < C{oi,m 2 ). (35) 

Proof. By virtue of (154]) 

IIVGH, < Ca||/i0||, + a||H| g + II div(K(g)gv (g> v) \\ q + //||AV x A|]g. (36) 
For q < 7 may we certainly write that 

at||/i<7|| g + ot||^u|| g < Ca;||f ||i 5 2 < Ca 3 ^ 2 , 



part are given 
(33) 
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by the continuity equation, the estimates (1141) and we get 



|| div(K(g)gv <g) v )\\ q < \\K(g)gv ■ Vv\\ q + a\\hK(g)v \\ q + a||^|| g 

< Cm 2 \\Vv\\ 2 q + Ca 3/2 , 

thus, the estimate (1241) of ||Vu|| g for 7 > q > 2 leads to 

|| dw(K[g)^v <g> < C(a 3/2 + a 2 + a 3 + a Q/q m\ +2{1 - 2/q) - 1 ) 
The last term in (1361) is bounded by the same constant, since 



||AVM||, < || Vu\\ q < a\\hg\\ q + a\\gv\\ q + || div(K(g)gv ® v)\\ q + 

+ C\\v ■ rlli-i/^+^an, 

where u is a weak solution to ( 12 5 p with a corresponding boundary condition after 
passing with e to 0, i.e. it satisfies 



— /jAuj = —arot(hg — gv) — rot div(K(g)gv ® v ) m 
w = \ 2x - w • r at <9fl 

□ 

For q such small that 7 > 1 + 2(1 — 2/0)7 we have then proved that 

Halloo < C{a)ml~\ (37) 
with 5 = 7(4/0 - 1) - 1 > and C(o;) = 6/g 

We will now apply the analogical decomposition for the approximative system (Q, i.e. 

v e = V0 e + V ± A e . 
Similarly as previously this leads to relation 

VG e = (2/i + u)A(j) e + P(g € ) 

= ahg — ag € v e — div(K(g e )g e v e <g) v e ) — eVg t Vv e + fxAV ± A e . (38) 

We are then able to prove that if e — > + the following lemma holds 

Lemma 7. G e — > G strongly in L 2 (f2). 

Proof. We will use the following fact: 

// V(G e — G) — ^ in L 2 , t/ien G e — G — >■ const in L 2 , 

and next we can show that at least for some subsequence e n — )• the constant is indeed 
equal zero 

\G e -G)= [ A(0-0 e )^O 
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since g£ = ^ = at dQ. 

on on 

This allows us to focus on showing the weak convergence, we have 
V(G e -G) = f iAV ± (A e - A) - a(g e v e - gv) 



- ( div(K(g e )g £ v e ® v e ) - div(K(g)gv <8> v)) - eVg £ Vv e . (39) 

The second term on the right hand side converges to weakly in L 2 owing to the strong 
convergence of v e — > v in L q for any < q < oo and by the boundedness of g e in L^. 
The last term converges to zero even strongly in L 2 . Now, by the continuity equation, 
the third term may be written in the form 



div(K(g e )g € v e <g> v e ) — div(K(g)gv <g) v ) = ahK(g e )v e — g e v e + eAg e v e 



+ agv — ahK(g)v + K(g e )g e v e ■ Vv e — K(g)gv ■ Vf , 

due to the argument explained above we need to justify the convergence only for two 
terms. Firstly note that eAg e v e converges to strongly in W / 2 _1 (fi). Secondly, since 
V(f e — v) — weakly in Lz(£i) we obtain the same information for K(g e )g e v e ■ Vf e — 
K(g)gv ■ Vv. 

In order to substantiate, taht the first term in fl39l) also tends to we observe that 

AV ± (A e - A) = V ± (cu t -cu), (40) 
and that the function u t — u satisfies the system of equations 



— jj,A{uj e — lo) = — avot{g e v — gv) — rot div(K(g € )g £ v £ ® v e — K(g)gv (g) v) 
—e rot(V^Vfe) in Q 

u e — u = ( 2x — — j (v e — v ) • r at dil. 

Repeating the same reasoning as in case of ( 126|) and above explications we can show, 
that V(o; e — to) consists of two parts. One of them converges to strongly in W 2 ~ 1 (Q) and 
the other converges weakly in L 2 (fi). Thus, by ([411]) . we get the same for AV ± (A e — A) 
and therefore the proof of lemma is complete. 

□ 

Provided with these information we can show the final argument for K(g) to be equal 
1 

Lemma 8. Let k > and let m satisfy 

||G||^ 7 < m < m 1 and \\G\\oo - ot(2/j, + v) > n > 

then we have 

lim |{x 6 O : Qe„( x ) > m }\ — 0- 
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Proof. 

The difference with respect to the Lemma 4.3 from |6j is that the rate of convergence 
here clearly must depend on a and thus we pass with e to when a is set. 
First observe that the assumptions of our lemma are satisfied. Indeed, as the difference 
m 2 — ||G||^ 7 increases with mi- Next, we introduce a function M(-) e C 1 (1R) given by 

{1 g < TO, 

g>m + l, 

G (0, 1) ge(m,m+l), 

where M'(g) < in (m, m + 1) and m + 1 < mi. 

We multiply the approximate continuity equation by M l (g e ) for some I G N and we 
observe 

a J M l (g e )(g-hK(g))dx + J M l (g e ) div(K(g)gv)dx = e J M l (g e )Ag dx 
n n n 

= -el J M'ig^M^ig^Vg^ 2 dx > 0. (41) 
h 

By integrating the second term on the left hand side by parts twice (the boundary 
terms disappear due to the definition of M(-)) one gets 

/ j J tM l - l {t)M'{t)dt J divv e dx 
n \o / 

>jj (hK(g £ )-g e ) dx+jj (g £ -hK(g e ))(l-M l (g £ )) dx. 

n n 

The first therm on the right hand side cancels due to the Theorem [3j We can replace 
divv e according to the definition of G e , then we have 

J ^jtM l -\t)M\t)dt^ (G e -P(g e )) dx 

a[2 ^ +v) J (g £ -hK(g e ))(l-M l (g £ )) dx. 



n \o 



< 



Since M'{t) is negative, supported in (m, m + 1) and m + 1 < mi — > m 2 the following 
inequality holds true 



m 



M l ~ 1 (t)M'(t)dt P{g £ ) dx 



n \o 



< m 2 



Qe 



M l - l {t)M'{t)dt 



\ Gt \ dx+ ^L±^l J \ ge -hK(g e )\(l-M l {g e )) dx. 
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After integration of the internal term we claim to conclusion that the above expression 
is different then only for a subset of Q, {g t > m}, thus 



m 
m 2 



J (l-M l (g e ))P(g e ) dx 

>m} 

J (l-M l (g e ))\G £ \ dx+ a ^ V) J \g £ -hK(g e )\(l-M l (g e )) dx. (42) 



{e t >m} 

< 



{g e >m} {g ( >m} 

Now for each 5 > we can find such sufficiently large number I e N, I = 1(5, e) that 

\\M l (g e )\\ L2({ee>m}) <5, (43) 

since M(g e ) is less then 1 for g e > m. This allows us to rewrite the inequality (1421) in 
the following form 

— — \{Qe >m}\ < —\\M l (g e )\\ L2i{s€>m}) \\P(g e )\\ L . 2{{ ^ >m}) 

ill 2 ill 2 

+ C(|n|)||G - G e || 2 + IIGIU \({g e > m}\ + a(2fi + v) \({g e > m}\ , 

where the term on the left is a consequence of the definition of P(-) and the limits of 
integration. By (j4"3"j) and the bound (125)1 we therefore may write 

f—-||G|| 00 -o ; (2 / / + \{g e >m}\ < —5c?l' 1 + C(\Q\)\\G - G e \\ 2 . 

Under our assumptions, the expression in the brackets is separated from and at least 
for a suitably chosen subsequence e n — > + Lemma [7] guarantees that 

lim \{g e > m}| < ^-^5a 3/2 . 

e n -*-o+ m 2 

As 5 may be arbitrary small and a = const, we truly have 

lim \{Qe n > m}\ = 0. 

□ 

This fact, as it was already mentioned before, completes justification that K(g) = 1 
a.e. in Q. 

The second problem to solve was to show that P(g) = P(g)- For this purpose we 
multiply the approximate continuity equation by the function In for 5 > and 
integrate over Q. Like in the proof of last lemma, we observe 

/ 777.2 f m 2 
In (g — h) dx + / In div(^)ax 

g e + 5 J g e + 5 

n n 

= e [ In -^^Ag dx = el [ dx > 0. (44) 

J + 8 J Qe + 5 

n n 
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Similarly as previously we integrate by parts, pass with 5 — > + , substitute G e from 
the definition and pass with e — > + to get 

J P{q)Q dx + (2/i + v)a J (g — h) In g dx < J Gg dx. (45) 
n n n 

From now on we will seek to reverse the sign of above inequality. We will use the fact 
that the limit continuity equation works with any smooth function up to the boundary. 
To indicate an appropriate one we first introduce the distribution: 

v ■ Vg = div(gv) — g divv. 

Then let us recall the following lemma (for the proof consult [7j). 

Lemma 9. Let Q G C ' 1 , v G Wg(Q), g G L p (0),l < p, q < oo, v ■ Vg G L S (Q), 
1/s = 1/p+l/q. Then there exists g n G C°°(Q) such that 

v ■ V g n —> v ■ Vg in L S (Q) and g n — >■ g in L p (fl). 

For such a g n one gets 



div(g n v)dx = J g n v ■ ndS = 0, 

n an 

thus passing with n — > oo our lemma provides that 

g divvdx = — J v ■ V gdx. 
n n 

Note that a function In — for 5 > is an admissible test function as it follows from 
the proof of Lemma |H] that < g n < m 2 , hence we get 



_ n 
Qn + 



n a 
We may now pass with n — > oo 

5 f gv ■ Vg 
g + S ' 



n n 

Next we also want to pass with 5 — > + , since j Q (g — h) In 5 dx = 0, the only difficult 
term is a J n h\n(g + 5), but it can be solved by the Lebesgue monotone convergence 
theorem, then we obtain 

a J hlng = a J glng — J v-Vg = a J ^ln^ + J g divv. 

n n n n n 
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Finally, recalling the definition of G one gets 



Gg dx = (2/i + v)a J (g — h) In g dx + J P(g)g dx. (46) 
q n 

The information contained in ( )45|) . (I46p together imply 

J P{q)Q dx + (2/i + z/)a ^ — h) In £ da; < (2/x + i>)a J (g — h) In g dx + J P{g)g dx. 
n n n n 

(47) 

The convexity of functions ghi.(g) and — h\n(g) ensure lower semicontinuity of the 
functional j n (g — h) ln(g) dx, in other words 



J (g — h)\n g dx < J (g — h)\ngdx. (4£ 



Therefore ( 14T|) reduces to 



P(g)g dx < / P(g)g dx. (49) 



By the 'standard arguments' we show that gg^ < g^ +1 which together with (141)1) yield 



e i e = e rfi. (50) 

Next, we may also show that g^ (x) < £> 7+1 (x) and g(x) < g^ 1 ^ for a. a. x G Q 

which easily imply that 

qI{x)=~(P(x). (51) 

Since L y (Q) is a uniformly convex Banach Space for 7 > 1, g e — £ weakly in L 7 (0) and 
II II 7 — ^ II f? II 7 we ma Y deduce, that g t — )■ £ strongly in L 7 (0). Thus in turn implies, 
that for some subsequence g t ^ g a.e. in f2 and the condition ||f? e ||L 00 (n) guarantees 
the uniform integrability of the sequence {£> en }^i which together with the Vitali's con- 
vergence theorem leads to the strong convergence of the approximate densities to the 
function g in L p (Q) for any 1 < p < 00. 



Remark 10. The density obtained in the above procedure is bounded by m as we could 
see in lemma\^ Now, by taking k sufficiently small and mi,m 2 sufficiently close to m 
the estimate for q — > 2 + with condition imposed on the assumptions of Lemma\E 
will imply that 

IMloo < C(a) 3 ^- 1 ), 

in particular, \\g\\oo < C(At)~ 3 . 
Theorem [T] is proved. 

□ 
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4 Passage with At — > 



+ 



In this section we wish to present the proof of Theorem El i.e. to demonstrate the 
passage with At — > + . The two previous section enable us to restrict attention to 
the case when the weak solution of the system (JTJ)- (j2J) exists, as provided by Theorem 
[TJ Our approach will be based on some estimates uniform with respect to the length 
of time interval At that we are going to gain here too. The task requires to work in 
the Bochner Spaces, but first let us introduce suitable notation: 

< t! X ' *J ~ t J^l / , , u *i j» w , if kAt < t < (k + 1) At. (52) 

This converts our original system into 

g[+ div(fw) = in Q, 

+ div(fw ®v)- yAv - (fi + u)V divv + Vn(g) = into, (f .^ 
v-n = atdtt, [ ' 

n- T(v, 7r)- r + /{>• r = at c?S7 

Moreover, the estimates (TTSJ) and ( I16p from the previous section now read: 

• g, g are bounded in L^O, T; L 7 (fi)) (54) 

• gv 2 , gv 2 are bounded in L^O, T; Li(O)) (55) 

• v, v are bounded in L 2 (0, T; tf 1 ^)) (56) 

• gv, 'gv are bounded in L^O, T; Ljj_ (fl)) U L 2 (0, T; L r (fl)) (57) 

7+1 

for 1 < r < 7 the last one holds as 

\\g k v% l/h+1) < ||e fc || 7 /2 ||f?V) 2 lli /2 and || Air < |b fe || 7 ||^||oo- £ , 
and all the bounds are independent of At. 

Our next aim will be to reconstruct the estimation for the norm of pressure Tt(g) = g 1 
in L q (Q x (0, T)) for some q > 1. Unfortunately, as we have seen in (j23p . such an 
estimate might not be true while q = 2, but it turns out to work for q — 1 + (I/7). To 
show this we test the momentum equation with a function $ of the form: 

$ fc = B{{g h ) - {g k }), into 
$ fc = atdVL 

From this testing we obtain the following identity: 

J (g k y +1 = J {g k V{g k }-J g k v k ®v k : V$ fc +^ J Vv k : V$ fc +(^+z/) ^ divv* div$ fc 
n n a n n 



+ /^(,v-A 1 )$ fc = E^ 

n i=1 
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Multiplying by At, summing over k = 1, . . . , M and employing our notation we get 



T T 

-7+1 




= / / (gy{g}-J J qv®v : V$+/i / / Wv : V$+(/i+z/) J J divv div$ 
on on on on on 

T 5 
+ / J ^ t (^-6(--^t)v(--At))^ = J2l i . (58) 
o n 1=1 

We go one with estimations for each of terms separately, 
(i) Since g is bounded in L^Li) and Loo(L 7 ) one gets 

T T 

1 




on o 
(ii) The Holder's inequality, (1561) and (1571) imply 

T T 



o n 



(in) Due to the properties of tha Bogovskii functional ||V$ fc || p < c(p, Q)\\g k \\ p , thus 



h + h = lij J Vv : V$ + {fx + v) J J divi)div<f>< j \\Vv\\ L2 \\V§\\ Ll+1 

o 

<CT^^\\ Q \\ L _ <+l{Li+1 y 



on on 



(iv) By the assumption that 7 > 2 we know that gv G £2(0, T; L2(0)) which is the 
special case of (157|) . hence by the continuity equation 





T 

i-(^-^(.-A^(--At))$ 
n 

T T 

+ J J -Lft - At)v(- - At)(<l(. - At) - $) 
on on 

T 

< sup [\g^\ + [ \\g(- - At)v(- - At)\\ L2{Q) \\g(t)v(t)\\ L2{n) 

0<t<T J J 



<c+ I U\\ljv\\l 2 ^_ 2) <c 
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All together leads to desired conclusion 

II^IC(l 7+1 ) < C (1 + T + T(- 1 )/ 2 ^)||^|| L7+l(i7+l) ) , 
in particular we have 

M 
k=l 

We are now in a position to validate that as At — > the following convergences 
hold: 

[g- g{- - At)}, [g- g] -> inL q (L 7 ) (60) 

for q G [1, oo) 

[gv — gv(- — At)], [gv — fw] — > in L q (L r ), (61) 
for {ge[l,oo), re[l,^]}U{ge[l,2), re[l, T )}, 

[fi)®t'-p'®ii]->fl m Li(L r ) n L q (Li), (62) 
for g G [1, oo) r G [1, 7). 

To see this it suffices to use the estimates (15^ [53| loTij 157)1 together with the observation 
derived from (j!5p . namely 

ll^-^(--At)||l 7(L7) <AtC, (63) 

moreover for the remaining term in (JT3j) we also have 

\\ g \v-v(.-At)\ 2 \\ Ll{Ll) <AtC. (64) 

From what has already been written we deduce that 

g, g — g weakly* in Loo(L 7 ), weakly in L 7+1 ((0,T) x f2), (65) 
v — ^ f weakly in L 2 (H 1 ). (66) 

Remark 11. Smce g g, v satisfy continuity equation f)53l) i . t/ius t/ie sequence of 
functions f(t) = (J n g4> dx) it) is bounded and equicontinuous in C[0,T] /or a// G 
C°°(0), • n = at <9f2. Therefore, the Arzela- Ascoli theorem, the density argument 
and the convergence established in yieid the following 

g, g g in C weak (L y ). (67) 

What is left is to show that we also have the corresponding convergence of the products 
gv, gv ® v. This can be done by repeated application of the following lemma. 

Lemma 12. Let g n , h n converge weakly to g, h respectively in L Pl (L P2 ), L qi (L q2 ) 
where 1 < p%,P2 < 00 and 

1 1 1 1 - 1 

Pi qi P2 <?2 

Let assume in addition that 
dg n 

-^j- is bounded in L\(W 1 m ) for some m > independent of n (68) 

\\h n - h n (- +i,t)\\ Lqi[Lq2) ^ as |f| 0, uniformly in n. (69) 
Then g n h n converges to gh in the sense of distributions on Q x (0, T) . 
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For the proof we refer the reader to [5]. 

For our case, since || is bounded in L 00 {W 2 ^^ +1 ^) and ^ is bounded in L 00 (W / 1 _1 ) + 
L 2 (H~ 1 ), the condition (I68p is satisfied for g n = g,^v and m = 1 respectively. Addi- 
tionally, we have that since h n = v is bounded in L 2 (H 1 ) the condition (l6T?j) also holds 
true. 

With this manner we see that gv converges weakly/ weakly* in Loo(L 27 /( 7+1 )) and in 
L2(L r ) for r G [1,7) to fw and that 'gv ^ v converges weakly in L\(L r ) D L q (Li), for 
g G [1, 00) r G [1, 7) to gv®v. Thus, the relations (161 j) and (I6"2"j) cause that we actually 
have 

£W — £W weakly in L q (L r ) (70) 
for {g G [1, 00), r G [1, ^]} U {q G [1, 2), r G [1, 7)}, 

fw ® — fw ® v weakly in Li(L r ) H L q (Li) , (71) 

for g G [1, 00) r G [1, 7). 

Having this we can pass to the (weak,weak*) limit as At — » + in the system (153|) 
everywhere expect in the term corresponding to the pressure: 

|f + div(gv) = in A, 

^ + div(gv tgiv)- fxAv - (/i + V divw + V7r(^) = in O, ^ 

n- T(t>, 7r)- t + /f • r = at dQ 

From now on we will be using the following denotation 

S(Vv) = n{Vv + V ± v) + v <iiv x vl, 
S(Vu) = /i(V{) + V-4) + z/ div^J. 

The proof of strong convergence of 7c(g k ) = (g k )' y in L^fi x (0,T)) is based on some 
properties of the double Riesz transform, defined on the whole IR 2 in the following way 

Kij - <IA A ),'<>,,. 

where the inverse Laplacian is identified through the Fourier transform J 7 and the 
inverse Fourier transform J 7-1 as 

(-A)-»=^ (|^F(tO). 

We will be using general results on such operators as continuity but also some facts 
concerning the commutators involving Riesz operators, being mostly a consequence of 
Div-Curl lemma [H] or that of Coifman-Mayer [T], [3J. The best overall reference here 
for both: auxiliary tools and the general idea of the proof is [3J. 

To take advantage of what we mentioned, there is a need to extended the system (153|) 
to the whole M 2 , as this is where the definition of the operator A" 1 makes sense. We 
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first observe that it can easily be done so for the continuity equation as gv ■ n = at 
dQ, hence 

^ + div(l^) = 0. (73) 
at 

For the momentum equation flo^ o we check that 

0(t,x)=i/>{t)((x)$, 4>= (V.A; 1 )^], 
^gC c °°((0,T)), (eCZ°(Q), 

is an admissible test function. This can be seen as a consequence of estimates (I54"t 155] 
[23 E3 E5]) and by the fact that the operator A" 1 gives rise to the spatial regularity 
to its range comparing to its argument of one. Particularly, later on we will take 
advantage of that for 7 > 2, the embedding W^(O) C C(Q) together with Remark [TT] 
imply 

(VsAj^lngl-KVA^lnf?] in C([0,T] xU). (74) 
Having disposed of this preliminary step, we can get the following integral identity 



where 



1 5 
J J VC(e 7 e + S(Vt)) : V.A-^lne]) cb dt = Y^Ii 

T 

h = J J tpC (^<9*0 + f?£ ® v : V a .Aj 1 V a ,[ln^|) da: dt, 



(75) 



Q 
T 



J J ^g J V x C • VxA-^ln^l dx dt, 



n 

T 



/ / ^S(Vv) : V.C ® VxAj^lne] dx dfc, 



n 

T 



y Jip{gv®v) : V x (® V^A^l^] dx dt, 



n 

T 



~ J J dtfPCoy ■ VxAx^lnol dx dt. 



n 

Analogically, if we test the limit momentum equation by the corresponding test function 
<p(t,x) = mC(x)<P, <\> = {VxA-^llng}, ^ E Cr((0,T)), C G C °°(^), (76) 
we get 

T 5 
y JipC(¥Q + ^v): V x A- 1 V x [l n g])dxdt = J2l i (77) 

n 4=1 
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where 



i 

h = J J tpC {ovdtt + qv®v: V^A^V^lf^]) dx dt, 



Q 
T 



o n 



y y ^S(Vw) : V,C ® V x A~ l [lne] dx dt, 



o n 

T 



J J ip{gv®v) : V X C® V.A^jl^] eft, 



o n 

r 



y y • v x A- i [i ng ] dx dt. 



o n 



The observation ( 171]) together with the consequences of lemma fT2l justify the conver- 
gences of the integrals I 2 , . . . , I5 from (175]) to their counterparts in (177]) . Thus we are 
left with the following identity 



lim / / ^C(g^~S(Vv) :K[l n o\) 

At— >0 




u 



— ]im^ (yQ v dt4> — QV ® v : 7?.[ln£?]J dxdt 



o n 

T 



J J 4>( {Q^Q - S(Vu) : H[lnQ\) dxdt 



o n 



y y (gvd t (f) — qv <g> v : 72.[lng]) dxdt. (78) 



o n 



By the continuity equation we obtain 

d t (j) = TZ[lugv], 

and the same we have for the test function in the approximate case, thus (178]) may be 
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rewritten as 

T 



lim Hj( {g'g - S(Vu) : K[lnQ\) dxdt 
At^oJ J 



o n 




— lim / / (qv1Z[1qqv] — gv <g) v : TZ[lng\) dxdt 

At— >0 11 



o n 

T 




V< {Q"<Q - S(Vu) : 7e[l n e]) dxcft 



o n 



~ J J^C (gvn[l a gv] -qv®v: TZ[l n g]) dxdt. (79) 
o n 

Now we will show that we actually have that the second terms on each of sides are 
equivalent as At goes to 0. For this purpose we will use the Feireisl's lemma [I] which 
is a consequence of the div-curl one. 



Lemma 13. Let 



V n — ^ V weakly in L p ( 
r n — r weakly in L q (I 



d2\ 



where 



1 1 1 

- + - = -< 1. 

p q s 



Then 

V n n(r n ) - r n 1Z(V 2 ) -± VTZ{r) - rK{V) weakly in L s ( m2 ' 



We will apply this lemma to r n = g(t,-), V n = gv(t,-) after extending them by 
on the rest of IR 2 and noticing that they satisfy assumptions of the lemma for p = 
27/(7 + 1), q — 7. Therefore we can take s = and thus, for a.a t £ [0, T) 

gvlZ(g)(t) — g1Z(gv)(t) — ^ gvlZ(g)(t) — gTZ(gv)(t) weakly in L S (Q) 

if we aditionally assume that 7 > 3. 

In view of this, the embedding L>2j_{Q) C W^i^l) and (J55|) we get that 

3+7 




hm^ / / tp(v (g7l[lngv] — gv1Z[lng]) dxdt 



n 



1 

//««(^[i„H-^[i n *^A, 



n 
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and the relations (ISTj) and (165 j) allow us to reduce ( I79p to 

T 

ton^ y (£ 7 £? - S(Vv) : TZ[l n g])dxdt 

T 

J y^C(F^-S(Vf) : K[l Q g\) dxdt. (80) 



o n 

T 



o n 



Now observe that by the fact that ( G C^°(Q) we may integrate by parts the second 
term on the left hand side and we will get 

T T 

f f 4>C S(Vv) : K[lng) dxdt = J JipK:[( S{Vv)) g dxdt 
on on 

T T 



on on 



(81) 



With the same manner we can transform the corresponding term in the limit on the 
right hand side of ( I80p . After passing with At to the limit in (18 ip we get 



^ J y V< S(V«) : n[l n g] dxdt 



o n 

T 



y y if>(2fM + v) Swvg dxdt + J j ' S(Vv)] - QK: [S(Vv)] , ) g dxdt 



on on 



32) 



where the precise form of last term on the right is a consequence of Div-Curl lemma. 
Therefore f[8"Dl reduces to 



i i 

j j 4< m - ^) «*«« = / / * m - e d*.«) dxA 

on on 
and since the choice of functions ^ and ( was arbitrary we have that: 



Q 1 Q — Q div x f = g"/g — g div x v. (83) 

Next, we take 5 > and multiply the discrete version of the continuity equation by 
ln(g k + 5). After integrating by parts over Q one get 

k Js-l\i„/Jc , x\ f „k„.k ^& 



n 

28 



By the Lebesgue monotone convergence theorem we can pass with 5 — > + and then 
integrate by parts once more to find 

■^ t J(e k - £ fc-1 ) Hs k ) + / div(«V = o- 

n n 

Recall that due to Theorem[T]we have J n g k = j Q g k ~ 1 , thus whereas x\n(x) is a convex 
function above equality may be changed into 

J [/ln(/) - /" 1 ln(/- 1 )] dx + J div{v k )g k < 0. (84) 



Summing from k = 1 to k = M, multiplying by At and using the notation ( I52p we 
transform (1841) to 



T 



J gln(g)(T) dx + J J g div^f) dxdt < J ^ln(^)(0) dx, 



n on 
thus after passing to the limit one get 

T 



J g\n(g)(T) dx + J J g dxv x v dxdt < J gla(g)(0) dx, (85) 
n on n 

For the limit momentum equation, we take advantage of the fact that it is satisfied in the 
whole space in sense of distributions, thus the solution is automatically a renormalised 
solution, i.e. it is allowed to multiply the equation by ln(^ + 5). Then we integrate 
over Q, pass to the limit as 5 goes to + and integrate with rescpect to time to get 

T 

glng(T)dx + J J g dxv x v dxdt — J g\ng(0) dx. (86) 
n on n 

By comparing the two results from (185]) and ( |86|) we get that 

T T 



g\n(g)(T) dx + J J g div x v dxdt < J g\ng(T)dx + J J g div^f dxdt. 
n on n on 

As in the proof of previous theorem, by the 'standard arguments' we show that gg^ < 
gi+ l which together with (183|) provide 

T T 




g di\ x v dxdt > J J g div^f dxdt. (87) 
on on 
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The two last information joined give the desired information, namely 



g In g = glng, 

and finally, by the convexity of function x In x, we obtain 

lim g = g a.e. in (0, T) x Q 

that completes the proof of Theorem [2j 

□ 
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